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Abstract - The classical method of power measurement of a 
system are iterative and bulky in nature. The new technique 
of measurement for bus voltage, bus current and power flow is 
a Phasor Measurement Unit. The classical technique and PMUs 
are combined with full weighted least square state estimator 
method of measurement will improves the accuracy of the 
measurement. In this paper, the method of combining Full 
weighted least square state estimation method and classical 
method incorporation with PMU for measurement of power 
will be investigated. Some cases are tested in view of accuracy 
and reliability by introducing of PMUs and their effect on 
variables like power flows are illustrated. The comparison of 
power obtained on each bus of IEEE 9 and IEEE 14 bus system 
will be discussed. 

Keywords - classical method; phasor measurement units 
(PMUs); state estimation; full weighted least square (WLS) 
state. 

I. Introduction 

The nature of electrical waves on grid will reflect the 
health of the system. The electrical wave is consisting of a 
complex number that represents both the magnitude and 
phase angle of the sine waves. The "Synchrophasor" is a 
device which measure nature of waves with respect to time 
[5] [ 10] . The PMU is a device that allows the measurement of 
voltage with multiple current options at each bus. The recently 
developed PMU [7] [8] [9] will help in deciding to stall power 
improvement devices at proper location for:- 

• Comprehensive planning 

• More accurate 

• Congestion tracking, 

• Advanced warning systems, 

• Information sharing, 

• Enhancement of System Integrity Protection 
Schemes (SIPS), 

• Quick restoration of Grids, 

• Effective grid operation, 

• Dynamically manage the grid. 

The classical method of measurement of power flow and 
voltage on a bus of power system are iterative and bulky in 
nature. The full weighted least square state technique is a 
nonlinear in nature, but with first order Taylor series it become 
a linear. Formerly conducted research will help in formulation 
of a relation between full weighted least square state technique 
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and classical technique. The natural approaches of parameter 
measurement will extravagance the PMU as computational 
burden on measurement. The optimal placement of PMU 
devices is a difficult task for power system state estimation is 
investigated, in the literature. This paper will reflect the effect 
of PMU on the accuracy of measurement on state estimation 
parameters. In case 1, the classical state estimation method 
without using any PMU. But in case P, the measurement of 
state estimation with PMU only is discussed. 

II. Wls State Estimation Method 
Consider the set of measurements given by the vector z:- 
Z= h(x) +e (1) 



Where: 
h T = 



[h, (x), h 2 (x), h, (x), h m (x)] 



(2) 



h (x) is the non-linear function relating measurement i to the 
state vector x 

x T = [x , x , x 3 xj is the state vector of system 

e T = [e, e, e, e lis the error measurement in state vector. 

L 1 2 3 m J 

The full weighted least square estimator [1] [3] [6] will 
minimize the following objective function: 



/(*)=£ 



(z,-W) 2 



[z-Kx)] T R- 1 [z-h(x)] 



(3) 



At minimum value of objective function, the first-order 
optimum condition to be satisfied. It can be expressed as 
follows: 

g (x) = ^^- = -H T (x)R- 1 [z-Hx)] = (4) 

ox 

The Taylor series of non-linear function g(x) can be 
expanded for the state vector x k by neglecting the higher 
order terms [2] will be as 

g(x) = g(x k ) + G(x k )(x -x k ) + = (5) 

The Gauss-Newton method is used to solve the above 
equation: 



[G(x k )]-> . g(x k ) 



(6) 



where, k is the iteration index, x k is the state vector at iteration 
k and G(x) is called the gain matrix and expressed as: 
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G(x) 



dg(x k ) 

dx 



H T (x k )R~ 1 H(x k ) 



(7) 



g{x k ) = -H T (x k )R- l [z - h(x k )] (8) 

Normally, the gain matrix is sparse matrix and decomposed 

into triangular factors. At each iteration k, the gain matrix are 

solved by using forward / backward substitutions, where 

Ax k+l = X k+l - X k and 

[G(x k )] Ax k+l = H T (x k )R~ l [z - h(x k )] = H T (x k ).RT l Az k (9) 
These iterations are going on till the maximum variable 

I A fc I 

difference satisfies the condition, ' Max \Ax \<S '. 



III. Classical Method 

The traditional/classical method of measurement is 
consider current as a relation of power flow with respect to 
bus voltages as 



(g^+b^iV'+V'-lViViCose^ 



v, 

The power injection at bus i can be expressed as, 

S^Pt + jQi 
G,=|v;E|v / |(G ( ,sm^-^cos^) 

.7=1 

The power flow from bus i to bus j are, 

S lJ =P iJ + JQ tJ 

ij=|Kr^+^)-|K||y;|^co4+^sin^) 

So the Jacobian H matrix will be as 



5P„ 



(10) 

(11) 

(12) 

(13) 

(14) 
(15) 
(16) 



H 



80 
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8P 

mi 


80 

sp^ 


30 
dO 

*~" in J 


80 

d Q flmv 


80 

mag 



8V 

dP flmr 


8V 
80 


8V 

d Q flmv 


8V 

mag 


8V 

8V 

mag 



8V 



(17) 



IV. Classical Method With Wls 

A PMU device will measure one voltage with multiple 
current phasors on a bus. Figure 1 shows a 4-bus system 
example which has single PMU at bus 1 . It consist of one 
voltage with three current phasor measurements, namely as 

V^e,, I, / 5 r I 2 Z8 2 and I 3 Z5 3 



Figure 1. Single PMU Measurement Model 

The transmission line normally defined as pie network 
due to their advantages on system constraints. If y is defined 
as the series admittance and y . as shunt admittance then 

J shunt 

the current measurements can be in rectangular coordinates 
as in Fig 2. 

Ii] = Cij -jDij 



AA/V^^"^ 



fc- -RHU*»- 



.-- KJ~Kk». 




Figure 2. Transmission Line in 'it' Model 

The expressions for current in real and imaginary 
component are as:- 

C, =| V, y„. | cos<$ +0J+\ Vj Y v | cos® +0,)-| V, Y„ 1 008(5, +^,)(18) 

D. =\V t Y*\sia® +0 sl )+\Vj Yj |sin(<5. +6j,)-|^|»n$ +6p(19) 
where, state vector is specified as: 

x = [| V, |Z0°, | V, \Z8 2 , \V 3 \Z8 3 \V N \Z8 N ] T (20) 

The Jacobian H matrix corresponding to their real and 
reactive parts is: 

-^ = |y„ |cos($+0j-|y e |cos(5,. + e tJ ) (21) 

dC 



dV 
dC. 



i=FJ«»(^+^) 
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= -\V t Y si |sin(<S, + 6J+\v t Y tj |sin(<S, + 9 U ) 
- = -\ Vj Y,\sm(S J + e,) 



(22) 



(23) 



(24) 
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dD, 



dV t 
dD 



dV 



f- = |y rf |sin(5 i + ^)-|F ( ,|sin(5 i + ^) 

|rJsin(<5, + 0,) 



dD, 



8S t 
dD 



OS 



VjYJcfXiSj + dg) 



(25) 



(26) 



(27) 



(28) 



The measurement vector z will be 

r ^ T n T ^ T It,! 7 " <?t 



T lT , 



z = t V > a; , v , gw • I y l • * ' c -j > d >j ] w 

Generally, measurements attain by PMUs are more 
accurate as compared to the traditional measurements. So 
that measurements performed with PMUs are projected to be 
more precise and accurate as estimated by classical methods. 

V. State Estimation With Pmus 
The state vector can be expressed as in rectangular 
coordinates. The voltage measurement (V = \V\Z. ) can be 

state as (V = E + jF), and the current measurement can be 
state as (I = C + jD). Where series admittance of the line as 
( g. + jb ) and shunt admittance of the transmission line as 
(g . ??+ jb . ). The flow of line current I can be expressed as: - 

i 9 =M -Vj)x(g ij +jb v )]+[VMgsi+JK)1 

= V,x [(g,j + jb tJ ) + (g si + jb si )]-VjX ( gij + jb tj ) 

The vector z is state as z = h(x) + e, (where x is a system state 

vector, h(x) is a linear equations, and e is an error vector). In 

rectangular coordinates: 

z = (H+jH m )(E + jF) + e (31) 

where, H = H+jH ,x = E + jF and z = A + jB. 

A and B are expressed by: 

A = HxE-H xF (32) 

B = H xE + HxF (33) 

m r v ' 

In matrix form, 



(30) 



H r 
-H. 



-H, 



(34) 



Then, the estimated value x = E + jF can be solved as:- 

AJc = {H T R l H) l H T R x kz = G l H T R x kz (35) 
If we define the linear matrix H„ as 

H r -H 
-H_ H.. 



H.,,. = 



then the eq. (35) can be written 



as: 
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X 



E 
F 



( H new R H ne«) H r, 



R 



(36) 



So that equation for x can be written in rectangular forms 
of z vector and H matrix. They are all in real numbers. In 
respect of the system, the PMU can deliver more precise 
information about system parameters. Some cases to be 
performed on classical measurement set with and without 
PMUs. The different cases simulations and analysis are as 
shown in Table I with some IEEE bus systems in the next 
section. 

Table I. Different cases in IEEE System 



Casej 


Measurements 


1 


Classical nith No PMUs 


2 


Classical with One IMU 


P 


Only PMUs 



VI. Simulation Results 

Some cases are tested for analyzing the system variables 
accuracy with or without PMU, with the help of MATLAB 
simulink software. The PMU locations in IEEE 9 Bus System 
and LEEE 14 Bus System at specific Bus Number are as shown 
in Table II. 

Table II. PMU Locations for Each IEEE System 



Typ e of Sys tern 


P1IU locations it Bus 


IEEE 9 System 


Bus 2.3 


IEE 14 Svstem 


Bus :.3 



The circuit diagram will be shown as in Figure 3 for IEEE 
9 bus system and Figure 4 for IEEE 14 bus system. 
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Figure 3. IEEE 9 Bus System 

In this segment, IEEE bus systems as IEEE 9 bus system 
[4] and LEEE 1 4 bus system [ 1 1 ] are tested with their respective 
cases to find out the consequences of the PMUs to the 
precision of the estimated variables. 
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Figure 4. IEEE 14 Bus System 

Assume n as the number of variables, m as the number of 
measurements and e as the ratio of the number of 
measurements per the number of variables. During the tests, 
maintained e as 1 .6. Table III has more detailed information 
about the measurement numbers for the tests [12]. 

Table III. Variable Type And Measurements 



Tjjies 
of 

SjslEm 


1 ar, 

(a) 


Mfia SIirEBlEat S 


Ratio 

E = 
(JH /») 


PowEr 


VoIcaeE 
Magnitude 


Total 
H 


IEEE 9 Bus 


17 


26 


1 


27 


lti 


IEEE 14 Bus 


27 


42 


1 


4? 


1.6 



A PMU has much smaller error deviations than classical 
measurements as 0.0000001 .The parameters measured are real 
power (flow & injected) and reactive power (flow & injected) 
measurements. The variation of parameters with or without 
PMU easily reflected in the Figure 5 - 12 as given below: 



The graph to be designed on the basis of Standard 
Deviation (S D) variation in each parameter's on each Bus 
which is planned for 20 cycle of operation of each IEEE 
system, it has standard deviation categorized in two 
categories i.e. Minimum & Maximum variation on actual 
parameters received at each Bus. The figure 5 & figure 7 will 
illustrate the variation in Standard Deviation at Minimum & 
Maximum of Real Power (P) with respect to Each Bus on 



■?MK 1 



P1IE 2 
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Figure 6. Graph P (SD) vs Bus Number 
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Figure 7. Graph P (SD) vs Bus Number 
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Figure 5 


Graph P (SD) vs Bus 
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Figure 8. Graph P (SD) vs Bus Number 
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Figure 9. Graph Q (SD) vs Bus Number 
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Figure 10. Graph Q (SD) vs Bus Number 



'QMAX 1 




-1.7SE-O0 ± 



BUS DUMBER 



2.O0E-03 



Figure 11. Graph Q (SD) vs Bus Number 

— ■ — QMAX 1 — ♦— QMA X F — *— q*I AX 2 

1.401-00 x 
1.2 0E-00 
1.00E-00 
p S.0 0E-01 



S.OOE-01 -- 
400E-01 ■■ 
2.00E-01 
0.00E-00 




A A A A A A 



T 7.00E-01 

6.00E-01 

-■5.O0E-01 

-■4.001-01 

3.00E-01!: 

2.00E-01 ' 

+ l.OOE-01 
O.OOE-00 
-1.00E-01 



12 J45f JB9 1011121314 
BUS NUMBER 

Figure 12. Graph Q (SD) vs Bus Number 

varying the number of PMU accretion from 'No PMU (Case 
1)' to 'Only PMU (Case P)' for IEEE 9 bus system. Similarly 
Figure 6 & Figure 8 will illustrate the variation in Standard 
Deviation at Minimum & Maximum of Real Power (P) with 
respect to Each Bus on varying the number of PMU accretion 
from 'No PMU (Case 1 )' to 'Only PMU (Case P)' for IEEE 14 
bus system. The Standard Deviation in Reactive Power (Q) 
will be illustrated in Figure 9 - 12 from 'No PMU (Case 1)' to 
'Only PMU (Case P)' for IEEE 9 bus system and IEEE 14 bus 
system. 

Table IV. Average S.D. In Estimated Variables 



Case 


\ s.Lut 


F 


Q 


9 Bus 


14 Bus 


.5 Em 


14 Bus 


1 


Mil 


-2.54E-4 


1.16E-5 


1.40E-2 


9.29E-3 


Max 


6.45E-4 


I.SBE-5 


-1.55E-1 


1.39E-2 


2 


Min 


-1.0BE-3 


-3.42E-2 


2 90E-2 


1.66E-1 


Max 


2.22E-3 


343E-2 


3.C7E-2 


B.04E-2 


P 


Min 


-137E-3 


-1.02E-1 


4.52E-2 


-2.B8E-2 


Max 


3.61E-3 


1.K5E-I 


4.77E-2 


2.10E-1 



The table IV illustrates that how the Standard Deviation 
(S D) enhanced at each stage as PMU increase from No PMU 
to Only PMUs. In IEEE 14 bus system, the S.D. of the 
estimated Power is approximately 1.88E-05 when there is no 
PMUs, but at Only PMUs it become 1 . 1 6E-0 1 . It means that 
the S.D. of 'No PMUs' is increased by adding PMUs in system. 
The interesting thing is that the standard deviation increasing 
as increasing PMU. Therefore, this result demonstrates that 
the effectively installing of PMUs will reduce the chances of 
error in measurement of estimated variables. 

Conclusion 

Now, the classical measurement method with the PMU is 
able to measure the voltage and current with their magnitude 
and phasors. The current measurement is implemented but 
the measurement set as in the rectangular form. The Jacobian 
matrixes are illustrated in detail for the measurements which 
include the elements of Equations. The state estimation in 
the linear formulation is investigated with PMU. All the 
variables and their respective measurements are improved as 
all in rectangular form, and then treated separately during the 
estimation process. Such linear formulation of the PMU data 
can produce the estimation result by a single calculation 
without performing the any bulky iteration as in classical 
methods. If only PMU data measurement set is exist in the 
real measurement world, and then there will be improvement 
in the computation time and accuracy estimation as compare 
to the linear formulation of the state estimation. 

The advantages of using PMU will advances the accuracy 
of the estimated variables. Some cases are tested while 
gradually increasing the number of PMUs which are added 
to the measurement set. With the help of advanced accuracy 
of PMU, it was seen that the estimated accuracy is also 
increases. One of the motivating thing is that the accuracy 
improves most effectively when the number of implemented 
PMUs are around '14 %' of the system buses. It is proved 
that the quality of the estimation is enhanced by adopting 
PMU data to the set of measurement. The PMU 
measurements will provide us improved accuracy and 
redundancy for the system. 

The study is carried out to establish a relationship among 
classical method, full weighted least square state estimation 
method and Phasor Measurement Units. Further it is visualized 
that the information of power variation on each bus of system 
are more precise and accurate as compared to the classical 
method. The linear formulation will suggest us a more specific 
and accurate information of power variation without doing 
any bulky iteration as performed in the classical methods. 
The use of PMUs will lead us to enhanced accuracy in results. 
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